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[J^ I Abstract. The present letter gives a rigorous way from quantum to classical random walks by 

introducing an independent random fluctuation and then taking expectations based on a path 
integral approach. 
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Sn 1 Introduction 

O ■ Recently the quantum walk (QW) is receiving significant attention, because it is a natural 
generalization of classical random walks to quantum systems, and because QWs could be 
useful in constructing efficient quantum algorithms. Here we focus on discrete-time QWs on 
p • the integers. QWs behave quite differently from classical walks. For the classical walk, the 
probability distribution is a binomial distribution with a variance that grows linearly with 
time. The probability distribution in a QW, by contrast, has a complicated and oscillatory 



^ ! form with a variance that grows quadratically with time. For recent comprehensive reviews, 
^^' see Kempe [T], Tregenna et al. [2], Ambainis jS]- 

A unitary matrix corresponding to the evolution of the QW is usually deterministic and 
^ ' independent of the time step. In this letter we study a random unitary matrix case motivated 
by numerical results of Ribeiro et al. jl] and Mackay et al. jnj ■ We obtain a classical random 
walk from the QW by introducing an independent random fluctuation at each time step 
and performing an ensemble average in a rigorous way based on a path integral approach. 
The time evolution of the QW is given by the following infinite random unitary matrices 
{Un:n = l,2,...}: 



Un 



K 



d. 



n 
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where an,bn,Cn,dn G C and C is a set of complex numbers. The subscript n indicates the 
time step. The unitarity of f/„ gives 



|2 I I |2 
^n\ T^ I '^n I 



where z is the complex conjugate of 2 G C and A„ = det Un = a^dn — &„c„ with | A„| = 1. 
Put Wn = {cLn, K, Cn, dn) ■ Let {wn : n = 1, 2, . . .} bc independent and identically distributed 
(or i.i.d. for short) on some space, (for example, [0,27r)) with 



E(aici) = 0. 



1/2, 



Remark that Eq.dlig) implies E{\ci\'^) = E{\di\'^) = 1/2, and Eq.dHSl) gives E{bidi] 
using Eq. p.ll) . The set of initial qubit states for the QW is given by 



(1.2) 
(1.3) 

Oby 



$ = {y. = *[«,/?] gC^ 



|aP + 



1}. 



where t means the transposed operator. Moreover we assume that {wn '■ n = 1, 2, . . .} and 
{a, 13} are independent. We call the above process disordered QW in this letter. Let R be 
the set of real numbers. Then we consider the following two cases: 

Case L a„, 6„, c„, rf^ G R (n = 1, 2, . . .) and E{aP + a/5) = 0. 

Case II: ^dap) = 1/2 and E{ap) = 0. 

We should note that Case I corresponds to an example given by Ribeiro et al. [^ and 
Case II corresponds to an example given by Mackay et al. p^, respectively. The numerical 
simulations by these two groups suggest that the probability distribution of the disordered 
QW converges to a binomial distribution by averaging over many trials. So the main purpose 
of the letter is to prove the above numerical results by using a path integral approach, which 
has been used in JHIIZIIHI, for example. In fact, our theorem (see Theorem 1) shows that the 
expectation of the probability distribution for the disordered QW becomes the probability 
distribution of a classical symmetric random walk. However our result does not treat a 
crossover from quantum to classical walks. Concerning the crossover, see jHl CHI HH 1121 UHl 
milTn], for specific examples. 

The rest of the letter is organized as follows. Section 2 gives the definition of the disor- 
dered QW. In Section 3, we prove Theorem 1. Section 4 is devoted to examples given by 



2 Definition of Disordered QW 

First we divide U„ into two matrices: 



an 



bn 




5 t^n 


" 

Cn 


" 

dn _ 



with Un = Pn + Qn- The important point is that P„ (resp. Qn) represents that the particle 
moves to the left (resp. right) at each time step n. By using P„ and (5„, we define the 



dynamics of the disordered QW on the hne. To do so, we define the (4A^ + 2) x (4A^ + 2) 
matrix f/„ by 



f/. 











Qn Pn 

g„ p„ 



. . . Qn 


... 



... g„ p„ 

Qn Pn 

Pn Qn 

where = O2 is the 2x2 zero matrix. Note that P„ and Qn satisfy 



„<^n 



Oo 



where * means the adjoint operator and I2 is the 2x2 unit matrix. The above relations 
imply that f/„ becomes also a unitary matrix. Let 



,(n) 



*rM 



,("), 






eC 



be a two component vector of amplitudes of the particle being at site k and at time n with 
the chirality being left (upper component) and right (lower component), and let 

be the state at time n. An initial state for the system is given by 



N 



2Ar+l 



v[/(o)(^)=*[0,...,0,<^,0,..::0]gC 

where = *[0,0] and ip = ^[a,P]. The following equation defines the time evolution of the 
disordered QW: for -A^ < k < N and I <n < N, 



For an initial state (/? = *[0, . . . , 0, (^, 0, . . . , 0], we have 

Af-l N-1 

U{^ = *[0~^, Pi(^, 0, Qi^, OT^, 

N-2 N-2 

U2UiTp = *[0r^, P2P1V, 0, iP2Qi + Q2Pi)^, 0, QsQi^, oT^- 

Let X^ denote the disordered QW at time n starting from ip E ^. By using Wf^\(p), the 
probability of {X^ = k} is defined by 



An) 



(n), 



,(n) 



pix- = k) = iiM^r^r = i^i:iMi^+ i^^:iMi 



Remark that {X!^ = k} is an event generated by {wi : i = 1,2, ... ,n} and {a,j3}. The 
unitarity of Um(jn = 1,2, ... ,n) ensures 



Y, p{x^ = k) = \\UnUn-i ■ ■ ■ UM? = w-ipf 



i«r+i/9r = i, 



k=—n 



for any 1 < n < A^. That is, the amphtude always defines a probabihty distribution for the 
location. 



3 Main Result 

We now give a combinatorial expression of the probability distribution of the disordered 
QW. For fixed / and m with I + m = n and —l + m = k,we introduce 



,{l,m) 



We should note that 



E 

L,mj>0:h-\ l-ln=l,mi-\ \-mn=m 



An) 



jln/^rn„ pln-i/^rn 



pin Qrrin pln-i Qrrin-i _ _ _ pliQrni 



hn™.! 



n^M = 2„(/,mV. 



For example, in the case of P{Xf = 0), we have 

S4(2,2) = P4Q3Q2P1 + Q4P3P2Q1 

+P4P3Q2Q1 + P4Q3P2Q1 + Q4P3Q2P1 + Q4Q3P2P1. 

Here we introduce the useful random matrices to compute 'Bn{l,'m): 

Rn 

In general, we obtain the following table of products of the matrices Pj,Qj,Rj and Sj 
(j = 1, 2, . . .): for any m,n > 1, 



Cn 




dn 




) »-'n — 
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■ 

hn _ 





Pn 


Qn 


Rn 


^n 


Pm 


O'm-t^n 


bmRn 


O-m-ttn 


h P 


Vm 


Cm^n 


(^m^n 


^m^n 


^m'^n 


Rrn 


Cm^n 


(^m^n 


Cni^n 


d P 


^m 


O^ra^n 


hmQn 


(^m^n 


'^m^n 



where PmQn = bmRn, for example. From this table, we obtain 
54(2,2) = 64(i3C2Pi + 040362(51 

+ (0463(^2 + &4C3&2)-Rl + (C4&3C2 + d4Csa2)Si. 

We should note that Pi,Qi,Ri and Si form an orthonormal basis of the vector space 
of complex 2x2 matrices with respect to the trace inner product {A\B) = tr{A*B). So 
'En{l,m) has the following form: 



S„(/, m) = pn{l, m)Pi + g„(/, m)Qi + r„(/, m)Ri + s„(/, m)Si. 



(3.1) 



In general, explicit forms oi Pn{l,fn),qn{l,m),rn{l,m) and Sn{l,Tn) are complicated. From 
Eq. ()3.1|) . we have 



S„(/,m) 



Pn{l, m)ai + r„(/, m)ci p„(/, m)6i + r„(/, m)(ii 
g„(/, m)ci + s„(/, m)ai g„(/, m)(ii + s„(/, m)6i 



Remark that p„(/,m), g„(Z, m), r„(/, m), s„(/, m) depend only on {wi : z = 2, 3, ...n}, so 
they are independent of wi, where Wi = {ai,bi,Ci,di). Moreover it should be noted that 
{wi : i = 1,2, .. .} and {a, P} are independent. Therefore we obtain 



P{X^ = k) 



||S„(/,m)v9||^ 

\\pn{l,m)ai +r„(/,m)cip + \qn{l,m)ci + s„(/,m)ai|^} \a\ 

+ ||p„(/,m)6i +r„(/,m)(iip + |g„(/,m)(ii + s„(/,m)6ip| | 
+ {(p„(/,m)ai +r„(/,m)ci)(p„(/,m)6i +r„(/,m)(ii) 
+ (g„(/,m)ci + Snil,m)ai){qnil,m)di + s„(/,m)6i)| a/3 
+ {(p„(/,m)ai +r„(/,m)ci)(p„(/,m)6i + r„(/,m)(ii) 
+ (g„(/,m)ci + Sn{l,m)ai){qn{l,m)di + s„(/,m)6i)| a/3 
Ci|a|2 + C2|/3|2 + Csa/? + C4a;5. 



Then we see that 



E{Ci\a\^) = -E[\p^{l,m)\^ + \sn{l,m)\^ + \qn{l,m)\^+\rn{l,m)\^]E{\a\^), 

since {p„(/,m), g„(/,m),r„(/,m), s„(/,m)} and Wi are independent, {wi : i = 1,2, . . .} and 
{a,/3} are independent, Eqs. ()1.2|) and ()1.3|) . Similarly, we have 

i?(C2|/3p) = ^E[\p^{l,m)\' + Kil,m)\' + |g„(/,m)|2 + |r„(/,m)p]i?(|/?|2). 

Case I: First we note that C4 = C3. So a„, 6^, c^, dn eH gives C3 = C4. Then the condition 
a/? + aP = implies Csa/^ + C4a/3 = 0. Therefore we obtain 

E{P{X^ = k)) 

= E(||S.(/,mV|p) 

= -E[\pn{l,m)\^ + \sn{l,m)\^ + |g„(/,m)|2 + |r„(/,m)|2] E(|a|2) 

+iE[K(Z,m)|2 + |s„(/,m)|2 + |g„(/,m)|2 + |r„(/,m)p]E(|/5p) 

= l{E(|p„(/,m)n+E(|g„(/,m)r) + E(M/,m)p) + E(|.„(/,m)p)}, 

since |ap + |/5p = 1. 

Case II: E{af3) = E{af3) = gives the above same conclusion in a similar fashion. 



To get a feel for the following main result, we will look at a concrete case of n = 4 with 
I = m = 2. 

E{P{Xt = 0)) 

= E(||S4(2,2V|p) 

= \ {E(b4(2, 2)n + E{\q,{2, 2)1^) + E{\n{2, 2)\') + E{\s,{2, 2)^)} 
= 1 {E{\h,\')E{\d^\^)E{\c,\^) + E{\c,\^)E{\a,\^)E{\h\') 

+E{\a,r}E{M^)E{\dp + i?(|64ni5|ic3ni5_a&2n 

+E{aih)E{b3cE)E{d2h) + E{a^h)E{hc3)E(d^b2) 
+E{\c,\')E{\h\')E{\c2\') + E{\d,\')E{\cs\')E{\a2\') 
+E{cd?jE{hci)E{c2a^) + E(clrf4)^fe)^(c^a2)} 

_ 1 1 2 2 _ 1 /4\ 
T6^T6^T6'^T6~24l^2J' 

since E{\ai\'^) = ^(|6ip) = ^(Icip) = E{\di\'^) = 1/2, and E{aici) = E{bid[) = 0. This 
result corresponds to P{Y^ = 0) = L)/2"^, where Y° denotes the classical symmetric random 
walk starting from the origin. We generalize the above example as follows: 

THEOREM 1 We assume that a disordered QW starting from ip satisfies Case I or Case 
II. For n = 0,1,2, ... , and k = —n, — (n — 1), . . . , n — 1, n, with n + k is even, we have 

Proof. For n = 0, 1, 2, ... , and k = —n, — (n — 1), . . . , n — 1, n, with n + k is even, let 
I = [n — k)/2, m = [n + k)/2 and M = [(n+k)/2)' ^^ general, as in the case of n = 4 with 
/ = m = 2, we have 

P{X^ = k) = E{\\E^{l,mM') 

-r M 

where u\ E {cj, bi, Ci, di} for any i = 2,3, . . . ,n, j = 1,2, . . . , M and Wk = {at, bk, Ck, dk) for 
any A; = 1,2,..., n. Then ^daip) = E(|6ip) = E{\ci\'^) = E{\di\'^) = 1/2 gives 

E{P{Xf^ = A;)) = ^ + E{R{wi, W2,..., Wn)). 

To get the desired conclusion, it suffices to show that E{R{wi,W2, . . . , w„)) = 0. When there 
are more than two terms in Pn{l,m) or Sn{l,m), we have to consider two cases: there is a A; 
such that 

• ■ ■ Pk+lPkPk-l ■ ■ ■ Pi, ■ ■ ■ Qk+lPkPk-1 ' ' ' -Pi, 



For example, 84(2, 2) = ^40302 + C463C2 case is k = 1, since 

Q4Q3P2P1, QaPzQ2Pi- 
Then the corresponding Pn{l,Tn) or Sn{l,iTL) are given by 



• ■ ■ dk+ldk^k-l ■ ■ ■ 0-2, 



■ ■ ■ Cfc+iOfcafc-i ■ ■ ■ a2- 



Therefore we have 



£'[(••• afc+iOfcafc-i ■ ■ ■ a2)(- ■ ■ 0^+10^0^-1 • ■ ■ 02)] 



E{- ■ ■) E{-a^Ck+i) E{\ak\^) E{\ak- 



E{\a2\ 



0, 



since E{aiCi) = 0. In other cases, a similar argument holds. So the proof is complete. 

The above theorem implies that the expectation of the probability distribution for the 
disordered QW is nothing but the probability distribution of classical symmetric random 
walk. 



4 Examples 



The first example was introduced and studied by Ribeiro et al. [4]. This corresponds to 
Case I: 



t/„ 



cosiOn 



sm 



sinyp 
- cosf 



where {^„ : n = 1, 2, . . .} are i.i.d. on [0, 27r) with 

E(cos2(^i)) = E(sin2(^i)) = 1/2, E(cos(^i) sin(^i)) = 0. 

Here we give two concrete cases: 

(i) di is the uniform distribution on [0,27r), 

(ii) P(^i = = ^(^1 = 7r/2 + = 1/2 for some ^ G [0, tt). 
For an initial qubit state, we choose a non-random y? = *[q;,/3] G ^E' with |a| = \(3\ = 1/v^, 
and aP + aP = 0. 

The second example was given by Mackay et al. [5]. This corresponds to Case II: 



Ur, 



1 

V2 



J9n 



where {^„ : n = 1, 2, . . .} are i.i.d. on [0, 27r) with 

E(cos(^i)) = E(sin(^i)) = 0. 

Moreover an initial qubit state is chosen as (^ = *[a, /?] G ^ with i5(|ap) = 1/2 and E{al3) = 
0. For example, both 61 and 6^ are uniform distributions on [0, 27r), and they are independent 
each other. Let a = cos{9^) and P = cos(6'*). So the above conditions hold. 



Finally we discuss an example studied by Shapira et al. jT3]. Let Wn = {Xn,Yn, Zn}, 
where {Wn : n = 1, 2, . . .} is i.i.d., and {Xn, Yn, Zn} is also i.i.d., moreover X„ is the normal 
distribution with the mean and the variance a^ with a > 0. Their case is 



^"^71 



1 1 
1 -1 



xK 



where 



K 



1 
V2 



cos(i?„)+zZ„^ 

,sin(_Rn) 



-Yn + iXn)- 



Rn 



{Yn + 
COs{Rn) 



■Y \ sm{R„) 
''^■^> Rn 

■y Sill(fin) 



Rn 



and Rn = J X"^ + Y^ + Z^. Then a direct computation implies 






Ei\h 



6 + 3 



1/2, 



'l_4^2x -2.2 



fi{(r) 



(4.1) 



We see that Eq. (jHH) implies < fiiV^/2) = 0.0179 . . . < /i(a) < l/2(= lim^|o/^(c^)) for 
any a > 0. Note that the limit a | corresponds to the Hadamard walk case. So their 
example does not satisfy our condition Eq. (jl.Hj) . 
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